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We present results for the 66 spectrum obtained using an 0(Mi,v )-correct non-relativistic lattice QCD action. 
Propagators are evaluated on SESAM's three sets of dynamical gauge configurations generated with two flavours 
of Wilson fermions at /3 — 5.6. Compared to a quenched simulation at equivalent lattice spacing we find bet- 
ter agreement of our dynamical data with experimental results in the spin-independent sector but observe no 
unquenching effects in hyperfine-splittings. To pin down the systematic errors we have also compared quenched 
results in different "tadpole" schemes and used a lower order action. 
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1. INTRODUCTION 

Lattice NRQCD has been widely used in the 
past few years to calculate a variety of phe- 
nomenologically interesting quantities. The first 
step in the NRQCD program, naturally, is the cal- 
culation of the spectrum of the 66 system. In the 
simulation presented here we use SESAM's large 
sample of dynamical Wilson-fermion gauge con- 
figurations to study both the spin-independent as 
well as the spin-dependent spectrum of the T. 
Our strategy, in searching for sea-quark effects, 
will be to compare our final dynamical results to 
that of a quenched simulation at equivalent lat- 
tice spacing. With our three sea-quark masses we 
can also study the dependence of mass-splittings 
on the light sea-quark mass. Recent simulations 
have exposed the sensitivity of spin-independent 
splittings to the details of the action. Therefore, 
in the spirit of we have implemented the 

NRQCD action including spin-dependent correc- 
tions of 0(MbV & ) and we tadpole improve using 
the mean link calculated in Landau gauge. With 
these ingredients we hope to clarify the effect of 
unquenching in the spin-dependent splittings. In 
addition, using the quenched configurations we 
investigate: (i) the effect of changing tadpole pre- 
scriptions; (ii) the effect of using an 0{M) } v i ) cor- 
rect action compared to the 0(Mi,v 6 ) corrected 
one. 
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2. SIMULATION DETAILS 

2.1. Action 

The non-relativistic (Euclidean) lattice Hamil- 
tonian to 0(Mf ) v 6 ) consists of S: the kinetic en- 
ergy operator, 

which is of order M^v 2 ; relativistic corrections of 
order M^u 4 , 

6HW = -ci ( A [ 2) } + c 2 -^L- ( A ■ E - E ■ A) 
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and spin-sensitive corrections of order Mf,v e 
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The E and B fields are represented by the stan- 
dard clover term or, where appropriate, by an im- 
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Table 1 

Simulation details 



p dyn = 5.6, n f = 2, 16^~x 32 



0.156 0.1570 0.1575 



# configurations 200 



200 



200 



measurements 



800 



800 



800 



/^quenched = 6.0, Uf = 0, 16 3 X 32 

# configurations: 200 
measurements: 800 

proved version to remove O (a 2 Mbi> 4 ) discretiza- 
tion errors. Following M the quark Greens func- 
tions are calculated from the evolution equation 



G(t+1) = 1- 
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with n = 2 |. We rely on tapole improvement 
choosing the mean link in Landau gauge as our 
improvement scheme. The coefficients Ci are then 
set to their tree level value of 1. 

2.2. Lattice Parameters 

The lattice parameters we have used are dis- 
played in table [l| Details concerning the gen- 
eration of our dynamical configurations and is- 
sues surrounding autocorrelations are discussed 
in rcf. ||. We exploit configurations more than 
once by starting the propagator evolution both 
at different spatial source points and on different 
timeslices. A binning procedure confirms that our 
4 measurements per configuration are indeed in- 
dependent. Throughout, we fix the bare heavy 
quark mass to a value aM^ = 1.7. 

3. SPECTROSCOPY 

For details concerning the choice of operators, 
smearing techniques and fitting procedures we re- 
fer to ||. We extrapolate energy level splittings 
linearly in the quark mass: 



Am = Amo + c 



(5) 



where m q denotes the bare light quark masses (see 
figure [j] for an example). The relevant scale we 



0.32 



0.3 
0.28 
0.26 
0.24 
0.22 

0.2 - 
0.18 - 
0.16 — 



-4-4 + 



CG( 3 P)-rS, 

0.01 0.02 0.03 0.04 0.05 0.06 0.07 



Figure 1. Extrapolation of spin independent 
splittings to m s /3.The triangular symbol denotes 
the extrapolated value. 





Average a 1 [GeV] 


i?SP 


n f = 0,p = 6.0 


2.47(10) 


1.43(9) 


k = 0.1560 


2.18(12) 


1.39(13) 


k = 0.1570 


2.32(10) 


1.33(12) 


k= 0.1575 


2.45(12) 


1.24(13) 


m s /3 


2.48(14) 


1.23(11) 



Table 2 

We use the average lattice spacing from the 
2 3 Sx - 1 3 S! and 1 3 P - l 3 ^! to convert our re- 
sults to physical units. Rsp is to be compared to 
the experimental value of 1.28. 



pick is given by m " + "^ +ms w 2i. The value of 
the strange quark mass is taken from our recent 
light quark mass calculation Q] and is such that 
Zjjf- ~ 0.0156 lies close to our lightest sea quark 
mass. 

Figure ^| and ^ summarize our results for 2 
and flavours. Unquenching effects are clearly 
visible in the spin-independent part of the bb 
spectrum. This is also evident from table |^: 
R SP = (2 3 S'i - l 8 £i)/(l 3 P - l 3 ^) disagrees 
with the experimental value i?sp = 1-28 in the 
quenched case but coincides with it when dynam- 
ical fermions are included. 

We do not observe any significant impact of un- 
quenching on the hyperfine splittings. In particu- 
lar, the P-hyperfine splittings seem to be underes- 
timated for both rif = and n/ = 2. Clearly, this 
result needs to be corroborated with higher statis- 
tics. Errors on our hyperfine splittings encompass 
the statistical error as well as the uncertainty in 
the fitting range; the latter is essentially of the 



3 



GeV 
10.5- 

10.0- * * 



9.5- 



tion as adding 0(Mi,v 6 ) corrections. 



3S 



2P 



2S IP 



IS 



i S 1 



'Pi 



Figure 2. The bb spectrum:radial and angular 
momentum splittings with the T-level set to its 
physical value. Open circles : n/ = 0,/3 = 6.0; 
filled circles : rif — 2,m q — m s /3; lines: experi- 
ment. 
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Figure 3. Fine structure: here the zero of energy 
is set to the T-level in the left part and to the spin 
averaged triplet P -level in the right part. Labels 
as in |2[ 

same size as the statistical error so that we can 
expect some improvement with higher statistics. 

We expect the inclusion of 0(M},v 6 ) terms to 
effect the hypcrfine splittings on the 10 % level, 
since v 2 « 0.1. Performing a relatively inexpen- 
sive 0(MbV 4 ) quenched simulation at (3 = 6.0 
(with 200 measurements) we find this naive ex- 
pectation to be well satisfied for the 3 Si — 1 So 
splitting as is shown in figure ^. However, using 
the plaquette prescription for uq (result from the 
NRQCD-Collab., also included in figure ||) shifts 
the value on the 10% level in the opposite direc- 




Figure 4. Hyperfine S-splitting obtained with an 
0(M b v 6 ) and C(M 6 w 4 ) correct NRQCD action. 
The higher order result is shifted significantly 
downward. Also shown: result with an 0(MbV A ) 
action but using the plaquette tadpole scheme. 
All results are for the quenched case with (3 = 6.0. 



4. SUMMARY AND OUTLOOK 

Overall, NRQCD has proven to work very well 
in the spin-independent sector, in particular in 
giving a precise determination of the lattice spac- 
ing a. It appears worthwhile to push the scale 
determinations to higher statistics using our con- 
figurations generated by T%L at k = 0.1575 on 
24 3 x 40 lattices (j3 = 5.6). This would enable 
a reliable extrapolation to 3 or 4 flavours. The 
dominant source of error in such a lattice scale 
determination is most likely due to the remaining 
lattice discretisation errors in the gauge configu- 
rations. 

Although higher statistics are highly desir- 
able, it seems unlikely that progress in the spin- 
dependent sector will come from this approach 
alone. It is worrying to find the spin-dependent 
corrections of 0(Mi,v 6 ) as large as 10 % even for 
the bb, decreasing the fine-splittings relative to 
the lower order estimates. A perturbative or non- 
perturbative calculation of the coefficients Cj is 
badly needed. 
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